log using "$Rep_smokelabor/2_analysis/output_logs/figure2.log", replace

	use "$Rep_smokelabor/1_build/regdata/county_day.dta", clear
	
	tsset countyfip rfrnc_dt
	
	reghdfe pm25 F(20(-1)1).hms_deep hms_deep L(1(1)20).hms_deep, a(county_doy st_yr) cluster(COUNTY10 rfrnc_dt) 
	summ pm25 if e(sample)
	
	** plot 
	import excel using "$Rep_smokelabor/2_analysis/output_figures/figure2.xlsx", clear sheet("pm25") first
	
	tw 	line l u eday, col(teal%70 teal%70) lp(dash dash) || ///
		connected b eday, col(blue) msize(0.8) ///
		xtitle("Days since smoke exposure") ytitle("Change in PM2.5 (ug/m3)") ///
		xline(0, lcol(black) lp(dot) lw(0.6)) ///
		yline(0, lcol(black) lp(dot) lw(0.6)) ///
		ylab(, nogrid) ///
		legend(off) ///
		graphregion(col(white))
	gr export "$Rep_smokelabor/2_analysis/output_figures/figure2.pdf", replace
	
log close 
